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Abstract 

We have investigated a contribution of the density of state (DOS) to the Mott metal-insulator 
transition (MIT). The MIT condition is given by the DOS, which usually requires the complicated 
fundamental equations. To simplify the process, we have derived a general solution of the MIT 
condition by an integral contains the DOS in the framework of coherent potential approximation. 
The formula is concise and presents the direct relationship between the MIT condition and the DOS. 
It is able to calculate the contribution quantitatively and to compare various DOSs analytically. 
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The electron correlation brings about various phenomena which are beyond the 

one-electron approximation. One of the most wide-ranging and historical issues is the Mott 
metal-insulator transition (MIT) 31. In order to discuss the MIT, for example, coherent 



potential approximation (CPA) 



6, 



11 



7 [, dynamical mean field theory (DMFT) 



0, y, Q, 



12] and their combination are used [l3j. For both of the methods, the single electron 



density of state (DOS) is an important factor for the MIT. The DOS is composed through 
all the eigenvalues of a model Hamiltonian for a single electron [l^. In this manner, the 
DOS is inseparably related to the Hamiltonian. Therefore, the condition for the MIT takes 
place varies by the DOS. For example, it can be qualitatively expected a material becomes 
a metallic state easily if it has more eigenstates near the band edge rather than around 
the band center. However, it is impossible to conclude immediately how the MIT condition 
depends on the DOS because the fundamental equations appear in CPA or DMFT are 
complicated ^, [l^. The numerical calculations, recent general approach, are also unable 
to give an answer for the question and the general dependence on the DOS. The numerical 
calculations are effective to treat the specified DOS. On the other hand, the analytical 
treatment is possible to handle the DOS in general form. 

Here, we have proven that the MIT condition is simply given by the DOS in the framework 
of CPA. We have derived the formula of the MIT condition for the symmetric DOS (e.g., 
it is valid for such as simple cubic and body centered cubic structures j^, [3]), in the case 
of paramagnetic state and zero temperature. The organization of this paper is as follows. 
First of all, we begin with a nondimensionalization of fundamental equations used in CPA. 
From the equations, we derive the necessary and sufficient conditions in order to constitute 
the energy band of electron. Using the conditions, we derive the MIT condition which is 
related to the DOS. 

To describe the Mott MIT using CPA, electronic Green function and the self energy for 
spin cr is determined self-consistently through these two equations 
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0,0, Q, 
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Here, u is energy of single electron, G{uj) is the single electron Green function and S(co') is 
the self energy. U is Hubbard U that represent repulsive potential between two electrons 
on the same site, n is average number of electron whose spin is —a. We treat n = 1/2 
for paramagnetic state. p'q{E) is normalized electronic DOS with no interaction between 
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them. After solving the simultaneous equations, we should consider for S +0. We begin 
our analytical treatment with the nondimensionalization of these fundamental equations 
([I]). Generally, p'^iE) >OatEo-D<E<Eo + D and p'q{E) = elsewhere. D 
is half width of the band width W, i.e., D = W/2. Introducing a dimensionless DOS 
Po{x) = Dp'q{E) and x = {E — Eq)/D^ J p'q{E)(\E = ^\ pQ{x)dx = 1 is obtained. For other 
variables, we define dimensionless variables as m = U/ a = {uj — Eq)/D — u/2, ai;{a) = 
ReS(u;)/D — u/2, b^ia) = ImS(co')/D, and e = 6/D. The dimensionless energy of single 
electron is expressed by a. ai;{a) and correspond to the real and imaginary part of the 
self energy, respectively (we omit the argument a). Also, we express dimensionless Green 
function as DG = Gr + iGj (Gr, G/ G -R). For simplicity, we omit the subscript e for Gr 
and Gj. Introducing = a — a^, we can rewrite the Green function as 
^ - x)po{x) /-i {b, - e)po{x) 

Thus, dimensionless fundamental equations are obtained as, 

Gnal + [1 - 2b,Gj] a, - (^b^ + ^^Gr = 0, G,a^ + 26,G^a, + 6, - (^6^ + ^ j = 0. (3) 

The hmiting values of a^, b^, and as e — +0 will be denoted by a, b, and ^ respectively 
if they exist. We define (a, b,^) as the solution of Eqs. ([3]). 

The electronic DOS with interaction between them is defined as p{a) = 
— l/vr lim£^+o > ojisl, Q. This definition requires 6 < 0. It is known if p{a) > 0, 

there is an electron which possesses energy a. This energy region forms an energy band. 
For example, if and only if m = 0, p{a) = po{a) is obtained and there is a single band. In 
order to discuss the Mott MIT, we need to figure out the energy band for u > 0. By the 
definition of p(a), the solution of Eqs. that satisfies p{a) > is confined to the next two 
cases: (i) 6 7^ and (ii) 6 = 0, |^| < 1. In the latter case, we obtain 

(^a'-^^Gn + a = 0, (^a'-^^G, = 0. (4) 

Gj ^ requires a = ±u/2 and obviously, Eqs. (jl]) fails for u > 0. Therefore, for all m > 
and a, p{a) > is equivalent to 6 7^ 0. If all the limiting values a, b, and ^ such that 6 7^ 
exist, Gr and Gj are also bounded and there are no indeterminate forms in Eqs. ([3]) in the 
case of e +0. In this situation, {a,b,C,) satisfies 

Gna^ + [1 - 2bGj] a - (b^ + ^ J Gr = 0, Gja^ + 2bGRa + b- (b^ + ^jGi = 0. (5) 



The Green function is 

We multiply Gi to the first of Eqs. ([5]) and Gr to the second of Eqs. ([5]). Subtracting the 
former from the latter, we obtain {2hG\ + 26Gj — G/)a = —hG^. Therefore, if (a, 6, ^) is the 
solution, it must satisfy either of below two sets of equations. The first set is given by 



-hGn 2 -4&(&G'| + hG] - Gi) [(1 - 2hGif + Ah'^G 
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(7) 



The second set is given by 



Gij = 0, 1 - IhGi = 0, + 6^ = — . (8) 

These two sets are necessary conditions for (a, 6, ^) is the solution of Eqs. ([5]). And obviously, 
if (a, 6,^) satisfies either Eqs. ([7]) or Eqs. ([8]), it is the solution of Eqs. ([5]). It means, these 
two sets are necessary and sufficient conditions for (a, 6,^) is the solution. Therefore, in 
order to constitute the energy band, we solve Eqs. ([7j) and Eqs. ([8]), instead of Eqs. ([5]). 

In the next, we consider the MIT in the case of p^^ix) = pq{—x). We define Uc as a critical 
value of u when the single energy band split into two regions. Namely, u < Uc represents 
a metallic state (the single band) and u > Uc represents an insulating state (the divided 
bands). We derive the expression of Uc analytically. 

To clarify our process of deciding Uc-, we determine the energy a that the splitting of 
the single band takes place. Because Pq{x) is even function, GR{b,C,) = —Gji{b,—^) and 
Gi{b,^) = Gi{b,—C,) are satisfied. As is obvious from Eqs. and (IHD, if (a, 6,^) is the 
solution, {—a,b, is also the solution for fixed u. Therefore, if p{a) > for a = a + ^, 
p(— a) = p{a) > is also valid for —a = —a — ^. This shows the energy band is symmetric 
with respect to a = 0. It makes a point that the band splitting takes place at a = 0. Thus, 
all we have to do is to find u that brings p{a = 0) > 0. Its supreme value is Uc- We search 
for all the possible u in Eqs. ([7]) and Eqs. (IHl). 

At first, Eqs. ([H]) fails for ^ 7^ 0. This is because for ^ 7^ and Gr = 0, we can calculate 

as 

1 _ 2bG, = 2 r 'f "f "^,T. .A - > »■ (9) 

Jo {x-^y + b^ {x + ^y + b^ ^ ' 

It means, Eqs. ([8]) can be satisfied only if ^ = 0. For ^ = 0, Gr = is automatically 
valid because Pq{x^ is even function. In this case, there is an unique 6 < that satisfies 
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1 — 2bGi{b,^ = 0) = and we can obtain a = a = -L'^u^ /A — W'. For a = 0, m = 2\h\ should 
be satisfied. As we show shortly after, this solution can be consolidated in the solution of 
Eqs. ([7]). Now, we turn on Eqs. ([7]). We rewrite a = a + ^ through the second of Eqs. ([7]) 
and obtain 

« = ^, o ^ X 



2hGl + 2hG] - a I Jo [{x - 0' + b^] [{x + 0' + b^] [iv - 0' + [(?/ + 0' + ' 

(10) 

It shows a = is equivalent to ^ = 0. With the fact that ^ = brings Gr = and 
1 — 2bGi 7^ 0, the last of Eqs. ([7]) becomes 



u 



b']- (11) 



\ \Gjib,^ = 0) 

Considering Gi{b,C, = 0) = 1/26 in Eq. (iTTl) . u = 2\b\ is obtained. This is the solution we 
derived from Eqs. ([8]). If 6 satisfies Eq. ( ITTl) for a given u, (a, b,^) = (0, b, 0) is the solution 
brings p{a = 0) > 0. On the contrary, for all 6 < 0, p{a = 0) > is obtained at the value 
of u which is given by the right side of Eq. (ITT]) . Therefore, we can formalize Uc as 



Uc = sup 

6<0 



A ^ I rl bpo{x)dx ^ I ' (l^) 
\ \J-1 ^2+^2 



This right side is monotonically decreasing function for 6 < 0. Since b —0 brings u — > +0, 

p{a = 0) > is obtained for u that increases continuously from 0. Consequently, taking 
b —CO, Uc is given by 

Uc = x^po{x)dx. (13) 

In conclusion, we have derived the Mott MIT condition with the DOS through the equa- 
tion, 

E- = . 2 [\^po{x)dx. (14) 



W \ Jo 

This equation is obtained under paramagnetic state, zero temperature and the DOS is even 
function in the framework of CPA. We have proven the dependence of the MIT condition on 
the DOS concisely and quantitatively. Equation (IT^ asserts that, if the integral becomes 
larger, the condition is satisfied with smaller W. The weight makes a larger contribution 
to the integral at |a;| ^ 1 than |x| ~ 0. It indicates a material metalizes easily if its 
eigenstates are collected near the band edge rather than around the band center. Also, the 
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integral of Eq. (fT^ is not more than 1. So in any cases, U <W is a necessary condition for 
the MIT. Through this proof, particularly in the key Eqs. iQ and ffTOl) . all the properties 
of po{x) we have used are just po{x) = po{~x) and po{x) > 0. Thus, to execute Eq. (IT^ . 



ta 



unctions 



14l |) and there is no 



16l . l20| . Once all the eigenvalues 



we can use the DOS poix) as its original form (sum of del1 
need to approximate the DOS by a continuous function 
for single electron is obtained, po{x) is composed and the MIT condition can be calculated 
analytically. 
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